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Abstract 

The action-angle variables for A^-particle Hamiltonian system with the Hamilto- 
nian H = J2n=o In s/i"^ (Pn/2) + In {p{xn - Xn-i) - p{xn + Xn-i)) , XN = xo, are 
constructed, and the system is solved in terms of the Riemann ^-functions. It is shown 
that this system describes pole dynamics of the elliptic solutions of 2D Toda lattice 
corresponding to spectral curves defined by the equation w"^ — P^{z)w + A^^ = 0, 
where P^{z) is an elliptic function with pole of order N at the point z = 0. 
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1 Introduction 



The main goal of this paper is to construct the action-angle variables for a finite dimensional 
Hamiltonian system of equations 

Xn = ixl-l){V{Xn,Xn+l) + V{Xn,Xn-l), Xn+N = Xn, (l-l) 

where 

V^u,v) = C(« - V) + an + V) - C(2.) = - i P> -")-P'e" + -") , ,1.2) 

2 p{u -V)- p{u + V) 

and to identify it as an elliptic analog of N -periodic Toda lattice. Here p{x) = p{x\2uj , 2uj') 
and ({x) = ({x\2u,2uj') are classical Weierstrass functions. 

Recently, finite-dimensional integrable soliton systems have attracted very special inter- 
est due to their unexpected relations to the theory of supersymmetric gauge models. The 
celebrated Seiberg-Witten ansatz ([|], 0) identifies moduli space of physically non-equivalent 
vacua of the model with moduli space of a certain family of algebraic curves. In [|], ^ it was 
shown that the family of curves corresponding to four-dimensional N = 2 supersymmetric 
SU{Nc) theory is defined by the equation 

Nc-l 

- wPnXE) + A'"^^ = 0, PnXE) = E''^+J2 (1.3) 

i=0 

In it was noted that this family can be identified with the family of spectral curves 
of Ai'c-periodic Toda lattice, and the Seiberg-Witten ansatz was linked with the Whitham 
perturbation theory of finite-gap solutions of soliton equations proposed in [§, 0] . Integrable 
systems related to various gauge models coupled with matter hypermultiplets in various 
representations were considered in p|-|3^, where more complete list of references can be 
found. 



In [^, ^ Ai'c-periodic spin chain related to XYZ model was proposed as soliton coun- 
terpart of iV = 1 supersymmetric SU{Nc) theory in six dimensions compactified in two di- 
rections, and coupled with Nf = 2Nc matter hypermultiplets. Spectral curves of A^c-periodic 
homogeneous XYZ spin chain have the form 

w'-wPliz)+Qf^^{z)=0, (1.4) 

where P^^{z) and Q^nS^) ^^re elliptic polynomials, i.e. elliptic functions with poles of order 
Nc and 2Nc at the point z = 0. Note, that ( p.. 41) is an elliptic deformation of the family 



of curves found in |3|] for four- dimensional N = 2 sypersymmetric SU{Nc) model coupled 
with matter hypermultiplets. 

A particular case of (|1.4|), when Qf^^^z) is a constant, Qfisf^{z) = A'^'^" can be seen as an 
elliptic deformation of (|1.3| ). The corresponding family of curves depends on A^^^ parameters 
which can be chosen as A and the coefficients Ut of the representation of P^{z) in the form: 

(_'[)N N-2 

P'n{z) = (^TZ^^f "V(^) + E u,dl-'p{z) + «o, (1.5) 
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An attempt to find a soliton system corresponding to the family of spectral curves defined 
by the equation 

w'' -wP^{z)+A'''' = 0, (1.6) 

led us to ( |lTD . After the system was found it turned out that, by itself, it is not new. Up to 
a change of variables g„ = p{xn), it coincides with one of the systems listed in where the 



classification of all Toda type chains which have Toda type symmetries was obtained. The 



new results obtained in this work are: isomorphism of with a pole system corresponding 
to elliptic solutions of 2D Toda lattice, the construction of action-angle variables, and explicit 
solution of the system in terms of the theta-functions. 



In [|T^, it was shown that a wide class of solutions of the Seiberg-Witten ansatz can 
be described in terms of a special foliation on the moduli space of curves with punctures. 
That allows one to consider such systems as reductions of 2D soliton equations. Following 
this approach, let us note, that (|L6|) defines an algebraic curve F as two-sheeted cover of 
the elliptic curve Fq with periods 2uj, 2uj'. Let P± be preimagies on F of 2 = 0. According 
to the construction of [Q, any algebraic curve with two punctures generates a family of 
algebro-geometric solutions of 2D Toda lattice 

(dl - 5L)^n = 4 (e^"+-^" - e^-'^"-) , (1.7) 

parameterized by points of the Jacobian J(F) of the curve. 

In the next section we show that algebro-geometric solutions (fn{x,t) corresponding to F 
defined by (|1.6|) are periodic in n up to the shift, ipn = fn+N + 2A^lnA, and have the form 

a{x - Xn+i{t) + a)a{x + Xn+iit) + a) 

ipn{x,t) = an{t) + In f.,:^r: ■ 

a[x — Xn{t) + a)a[x + Xn[t) + a) 

Substitution of (|1.8|) into ( |1.7|) leads to equations ( [L.l|) for x„(t). 

It section 3 we construct a new Lax representation for ( |1.1| ) and show that the spectral 
curve defined by the Lax operator has the form (|1.6|) . We prove also, that if Xn{t) is a 
solution of ( |1.1| ), then there exist functions a„(t) (unique up to the transformation a„(t) — > 
c^nit) + cit + C2, Ci = const), such that the functions (fn{x, t) of the form (|1.8|) satisfy (|1.7|) . 

The last section is devoted to the Hamiltonian theory of system (|1.1|). Equations ( |1.1| ) 
are generated by the Hamiltonian 

Af-l 

if = ^ In Sh~^ {Pn/2) + In {p{Xn - Xn-l) - p{Xn + Xn-l)) , (1-9) 
n=0 

and the canonical Poisson brackets {pm, Xn} = Snm- We would like to emphasize that though 
this Hamiltonian structure can be easily checked directly, it was found by the author using 
the algebro-geometric approach to Hamiltonian theory of the Lax equations proposed in 
T7| , and developed in The main advantage of this approach is that it allows us 



to find simultaneously the action-angle variables and a generating differential which defines 
low-energy effective prepotential. 

Note, that from the relation of system (|1 . 1| ) to 2D Toda lattice it is clear that degeneration 
of the elliptic curve Fq corresponds to a degeneration of this system to the Toda lattice. 
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It would be very interesting to consider this degeneration explicitly on the level of the 
Hamiltonian structure. We will consider this problem elsewhere. 



2 Elliptic solutions of 2D Toda lattice 



Algebro-geometric solutions of 2D Toda lattice were constructed in [^. Let F be a smooth 



genus g algebraic curve with fixed local coordinates z±{Q) in neighborhoods of two punctures 
P± G r, z±{P±) = 0. Then, for any set of g points 71, . . . , 7^, in general position there exists 
a unique function ipn{x,t,Q) such that: 

1". ilJn{x,t,Q), as a function of the variable Q G F, is meromorphic on F outside the 
punctures P± and has at most simple poles at the points 'js (if all of them are distinct); 

2". in the neighborhoods of the punctures the function ip^ has the form 

= 4^e(-±*)^"^ (Eef(x,t)4) , ^0+ = 1- (2.1) 

Uniqueness of ipn implies that it satisfies the following system of linear equations 

{dt + dr,)ilJn{x,t,Q) = 2iljn+i{x,t,Q) +Vn{x,t)ilJn{x,t,Q), (2.2) 

{dt- d.j:)lpn{x,t,Q) = 2Cnix,t)'llJn-lix,t), (2.3) 



where the coefficients are defined by the leading coefficient of expansion (|2.1| ) with the 
help of the formulae: 

Vn = {dt + d,)Mx,t), c„ = e'^"(^'*)-^"-i(^'*), ^n{x,t)=\n^^{x,t). (2.4) 

Compatibility of ( p.2|) and (|2^ ) implies that ipn{x,t) is a solution of 2D Toda lattice ([L^ ). 

The function ilJn{x,t,Q) is called the Baker-Akhiezer function and can be explicitly ex- 
pressed in terms of the Riemann theta-function associated with a matrix of 6-periods of 
holomorphic differentials on F. The corresponding formula for ip^ is as follows. 

Let us fix a basis of cycles a^, 6,, i = l,...,g,onr with the canonical matrix of intersec- 
tions: o aj = bi o bj = 0, Oj o bj = 6ij. The basis of normalized holomorphic differentials 
dQ'j{Q), j = 1, . . . , g, is defined by conditions dQ^ = 6ij. The 6-periods of these differen- 
tials define the Riemann matrix Bkj = §i,^ dVL^. The basic vectors Ck of and the vectors Bk-, 
which are columns of the matrix B, generate a lattice B m . The ^f-dimensional complex 
torus 

J{T)=C^/B, B = ^nkek + mkBk, nk,mkeZ, (2.5) 
is called the Jacobian variety of F. A vector with the coordinates 



AkiQ) = / dnt (2.6) 
Jp+ 



defines the Abel map A : F — > J{^)- 
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The Riemann matrix has a positive-definite imaginary part. The entire function of g 
variables z = {zi, . . . , Zg) 

e{z) = e{z\B) = J2 e2^*(^'"')+^*(^'"'"\ 

meZ9 

is called the Riemann theta-function. It has the following monodromy properties 

eiz + ek) = e{z), 9{z + Bk) = e-2--'=— (2.7) 

The function 9{A{Q) + Z) is a multi- valued function of Q. But according to ( |2.71) , the zeros 
of this function are well-defined. For Z in a general position the equation 6{A{Q) + Z) = 
has g zeros 71, ... ,7^. The vector Z and the divisor of these zeros are connected by the 
relation Z = — ^{is) + K^, where /C is the vector of Riemann constants. 

Let us introduce normalized Abelian differentials dVl^^^ and d^l^^^ of the second kind. 
They are holomorphic on F except at the punctures P±. In the neighborhoods of P± they 
have the form 

rffi(^) = d{z^^ + 0(1)), rffiW = d{±zg^ + 0(1)). 

Normalized means that they have zero a-periods. The vectors of 6-periods of these differen- 
tials are denoted by 27TiV and 27riW, i.e. the coordinates of the vectors V and W are equal 




Let dfi*^") be a normalized abelian differential of the third kind with simple poles at the 
punctures P± with residues =f1. From the Riemann bilinear relations it follows that the 
vector of its 6-periods satisfies the relation 

Uk = ^<f = A(P_) - A(P+). (2.9) 

Zm Jbk 

If we choose a branch of the Abelian integral OS"^' near P+ such that = -\nz+ + 0{z+), 
then near P_ it has the form 

^(n) =inz_+/o + 0(^_). 



Theorem 2.1 (ji^j) The Baker-Akhiezer function is equal to 

^ixtO)- (^i^MQ)+nU + xV + tW + Z)e{Z) , 

e{nU + xV + tW + Z)e{A{Q) + zf''^y'^^ +a;iZ +t\l ). [Z.iU) 

The function (fn{x,t) given by the formula 

, , , , ^((n + l)t/ + xF + tW^ + Z) 
X, t) = nio + In ^\ / 2.11 
0[nU + xV + tW + Z) 

is a solution of 2D Toda lattice. 
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For a generic set of the algebro-geo metric data the function ipn{x,t) given by ( p.llD is a 
quasi-periodic meromorphic function of all the variables {n, x, t). In the solutions of 2D 
Toda lattice which are elliptic in the discrete variable n were considered. It was found that 
dynamics of its poles coincides with the elliptic Ruijsenaars-Schneider system In this 
paper we consider solutions that are elliptic in the variable x and are periodic in n. 

The condition that is elliptic in one of the variables is equivalent to the property 
that the complex linear subspace in J{T) spanned by the corresponding directional vector 
is compact, i.e. it is an elliptic curve Tq. In the case of x- variable it means that the vectors 
2uaV, a = 1,2, belong to the lattice B defined by (|2.5|) : 

2uJaV = J2n'^ek + mtBk, nl,ml e Z. (2.12) 

k 

Here and below uoi = uj,uj2 = oj' are half-periods of the elliptic curve Fq. 

Theorem 2.2 Let T he a smooth curve defined by equation and let P± be preimages 

on F of the point z = eTq with local coordinates in their neighborhoods defined by the local 
coordinate z on Tq. Then the corresponding algebro- geometric solutions given by formula 
( \2.11\ ) satisfy the relation 

ipn+N{x,t) =<fn{x,t) + 2N In A, (2.13) 



and have the form (lj_S), i.e. 

a{x — Xn+iit) + a)a{x + Xn+iif) + a) 



(Pnix,t) = an{t) + In- 



a{x — Xnit) + a)a{x + x„(t) + a) 



The functions Xn{t) defined by this representation enjoy equations The functions an 

satisfy the relation 

4e""- = (l - Xl{t)) W{Xn, Xn+l)W{Xn, Xn-l), (2.14) 

where 

tnu.,,)- -^'" - ">;'"+ "' . (2,15) 

a[2u) 

Proof. The first statement of the theorem is a direct corollary of the uniqueness of the 
Baker- Akhiezer function. The projection Q = {w, z) & T ^ w defines w = w{Q) as a 
function on the curve. This function is holomorphic on F outside the puncture P+, where 
in has the pole of order N, w = z^^{l + 0{z)). At the point P it has zero of order A^, 
w = A^^2;^(l + 0{z)). Therefore, we have the equality 

ipn+Nix, t, Q) = w{Q)tpn{x, t, Q), (2.16) 

because the functions defined by its left- and right-hand sides have the same analytical 
properties. 

Let us consider the functions 

T«(z) = e2fW--2,..^ ^ ^2.17) 
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They are double-periodic and holomorphic on Tq except at 2; = 0. Again, comparison of 
analytical properties of the left- and right-hand sides proves the equality 



ipnix + 2uJa,t,Q) =Ta{z)ipn{x,t,Q), Q = {w , z) . (2.18) 

The function e"^" is defined as a ratio of the leading coefficients of an expansions of ipn on two 
sheets of F. Therefore, it does not change under the shifts x — ^ x + 2uja, and consequently, it 
is an elliptic function of the variable x. From ( |2.11| ) it follows that if we denote roots of the 
equation 9 {nil + xV + tW + Z) = in the fundamental domain of Fq by x;^ (t) , j = 1, . . . , D, 
then 

j^i a{x - xi{t)) 

Our next step is to show that e'^^ has only two poles and zeros in Fq. 

Lemma 2.1 The function 6{xV + ^) corresponding to a smooth algebraic curve F defined 
by as a function of the variable x is an elliptic theta-function of weight 2, i.e. it can 

be represented in the form 

e{xV + = r{0(T{x - x\OMx - x\0)- (2-20) 

Proof. Let us find the coefficients of expansion ( |2.12| ). The branching points zf of F over 
Fq are roots of the equations Pfl{z) = ±A^. In a generic case, when they are distinct, the 
curve F is smooth. The Riemann-Hurwitz formula 2g — 2 = u which connects genus g of 
branching cover of an elliptic curve with a number u of branching points, implies that F has 
genus + 1. We choose a^, hi cycles on it as follows: 

tti, i = 1, . . . , N — 1, are cycles around cuts between branching points zf ^z^ and qn and 
flAT+i are two preimages of a-cycle on Tq. (We assume that a and 6-cycles on Fq correspond 
to the periods 2uj and 2uj', respectively.) 

From the definition of the differential d^l^^^ it follows that 

= d (^C{z) - . (2.21) 
Therefore, the coordinates of the vector V defined by ( p.8| ) are equal to 

V, = 0, t = l,...,N-l, Vn = Vn+i = -(v'--^')=-i^- (2.22) 

m \ UJ J 2uj 

Comparing the vector of 6-periods of dVt^^'^ with the vector (0, . . . , 0, 2uj\ 2uj') of 6-periods of 
the differential dz, considered as a differential on F, we get 

= — I dz, i = l,...,N + 1. (2.23) 

hi 2uJUj' Jbi 

The a-periods of dz are equal (0, . . . , 0, 2uj, 2u). Therefore, 



dz = 2uj{dn% + dn^ 
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where dVL^ are normalized holomorphic differentials. From (|2.23| ) we finally obtain that 

2uj'V = -Bn - Bn+i, (2.24) 
where Bi are the vector of 6-periods of dTl^. The monodromy properties of 6'-function imply 

e{{x + 2uj)V + Z) = 6{xV + Z), 6{{x + 2uj')V + Z)=e^^''^e{xV + Z), (2.25) 

where 

l{x) = Txi (2a;(VAr + Viv+i) + Bn+i,n+i + Bn^n — Bn,n+i — Bn+i,n + 2Zn+i + 2Zn) 
Using (|2.22| ) we obtain 

dUx) = dx. (2.26) 

CO 

The number D of zeros of the function 6{xV + in the fundamental domain can be found by 
integrating of the logarithmic derivative of this function over the boundary of the domain. 
From (|2]2|) and (|]26D it follows that 

D = -1-<X d\n9(xV + Z) =2. (2.27) 

2m JdGo 

The equality (|2.2CI|) is proved. It implies that the index j in ( p.l9|) takes values j = 1,2. The 
sums of zeros and poles of an elliptic function are equal to each other (modulo periods of 
Fo). Hence, x{(t) can be represented in the form 

xl, = Xn{t) + a{t), xl{t) = -Xn{t) + a{t). (2.28) 



In order to complete a proof of (|1.8| ) we need only to show that a(t) does not depend on t. 

Let us substitute ( |2.19| ) into (|1.7D. A'priory the difference of the left- and right-hand sides 
of ( |1.7|) is an elliptic function of x with poles of degree 2 at the points x{(t) and Xn_^^i(t). 
Vanishing of the pole of degree 2 at implies that 

(i;y'-l = F:«), (2.29) 

where 

Uj (^{x - xi^^)a{x - xi_^] 



F'n{x)=r^ ' ^ "-^^ , r„ = -4e--— . (2.30) 

lljjti \X — Xr. 



Vanishing of the pole of degree 1 at implies 

(Ec«-^^n4-i) + c«-4-i)-2Ec«-^f.)| • (2.31) 

^ V i jj^i J 

Substitution of (|]28D in (ICTI) shows that 

F^{xl) = Fl{xl) = rn{t)W{Xn,Xn+l)W{Xn,Xn-l). 

Hence, we obtain the equahty {xD^ = {x'^Y, which implies that d = 0. Equalities ( |1.8|) and 
( p.l4| ) are proved. At the same time substitution of ( |2.28| ) into ( p.31[ ) gives equations (|l.lj). 
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3 Generating problem and Lax representation 



In this section we construct the Lax representation for (|1 . 1|) following an approach proposed 
in and developed in ^ ^ (see their summary in |^)). According to this approach 



pole dynamics can be obtained simultaneously with its Lax representation from a specific 
inverse problem for a linear operator with elliptic coefficients. 

In the most general form this inverse problem is to find linear operators with elliptic 
coefficients that have sufficiently enough double-Bloch solutions. A meromorphic function 
f{x) is called double-Bloch if it has the following monodromy properties: 

f{x + 2uo.) = BJix), a = 1,2. (3.1) 

The complex numbers Ba are called Block multipliers. (In other words, / is a meromorphic 
section of a vector bundle over the elliptic curve.) It turns out that existence of the double- 
Bloch solutions is so restrictive that only in exceptional cases such solutions do exist. 

The basis in the space of double-Bloch functions can be written in terms of the funda- 
mental function $(x, z) defined by the formula 

$(a;,.) = 4^e«^)^ (3.2) 

From the monodromy properties of the Weierstrass functions it follows that $, considered 
as a function of z, is double-periodic: $(x,z + 2ua) = $(x, z), though it is not elliptic in 
the classical sense due to essential singularity at 2; = for x ^ 0. As a function of x, the 
function $(x, 2;) is double-Bloch function, i.e. 

$(x + 2ua, z) = Ta{z)^{x, z), 



where Ta{z) are given by ( |2.17| ). In the fundamental domain of the lattice defined by 2000 



the function $(a;, z) has a unique pole at the point x = 0: 

$(x,z) = + 0(x). (3.3) 

Let f{x) be double-Bloch function with simple poles Xi in the fundamental domain and with 
Bloch multipliers Ba (such that at least one of them is not equal to 1). Then it can be 
represented in the form: 

N 

/(x)=^Q<l>(a;-x„^)e^^ (3.4) 

i=l 

where q is the residue of / at Xi and {z, k) are parameters such that Ba = Ta{z) exp{2uJak). 
Now we are in position to present the generating problem for ( |L1|) . 

Theorem 3.1 The equation 

{dt + d,)^n = 2^„+l + Vn{x, t)^n (3.5) 
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with an elliptic coefficient of the form 

2 

Vnix,t) = 7„(t) + J2 [hnit)C{x - xi{t)) - /i^+i(t)C(x - , (3.6) 

i=l 

where 

xl^{t) = Xnit) + CL, = —Xnit) + tt, tt = COUSt, (3.7) 

has two linear independent double-Bloch solutions with Bloch multipliers Ta{z) (for some z), 
i. e. solutions of the form 

vl/„(a;,t) = 5:<$(x-<(t),z) (3.8) 

i=l 

if and only if the functions Xn{t) satisfy equation 

// equation l \3.^ has two linear independent solutions of the form l \3.^ for some z, then 
they exist for all values of z. 

Proof. Let us substitute ( p.8| ) into ( |3.5| ). The both sides of the equation are double-Bloch 
functions with the same Bloch multipliers and with the pole of order 2 at x^, and the simple 
pole at xl^_^_i. They coincide iff the coefficients of their singular parts at these points are 
equal to each other. The equality of the coefficients at {x — x^)"^ implies 

h^ = in- 1- (3.9) 
The equality of residues at xl^_^_i is equivalent to the equation 

The equality of residues at is equivalent to the equation 
where 

K = 7n - E ^n+lC« - ^^n+l) + E ^'nC« ' <)■ (3-12) 

Equations ( p. 101) and ( p.ll|) are linear equations for c^. Their compatibility is just a system 
of the equations: 

dt{\^K+i)^{<+i - <) + (in+i - a:{)$'«+i - <) = (M;+i - Mi)^{x\^^ - + 

+ E - - ^n) - E - a^n)^''^ (a^n " 4), (3-13) 

ki^i k^j 

which can be written in the matrix form: 

dtL^ = M„+iL„ - L„M„, (3.14) 

where and M„ are matrices defined by the right-hand sides of ( p.lO| ) and ( |3.11| ). Equations 
( p.l4| ) are necessary and sufficient conditions for the existence of solutions of ( p.5|) which have 
the form ( |3.8| ). Therefore, the following statement completes a proof of the theorem. 
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Lemma 3.1 Let = [L^^it, z)) and Mn = {M^(t, z)) be defined by the formulae 

LH = 2-'hi_,M<+i = Ml Mi, = -xlz),i^ J, (3.15) 

where x\ = x„, x\ = — x„, h]^ = x\ — 1, and M^ is given by (^.1^ ) with 7„ such that 

{xl-l)a\2x^) \ 



In - 7n-l = dt In 



Then they satisfy equation (\3. 14 ) if and only if the functions solve equations 



(3.16) 



Note, that ( |3.16| ) defines 7„(t) up to a constant shift jnit) —>■ ln{t) + 9{t), which corre- 
sponds to the gauge transformation — > e^'ifn of equation (p.5|), and which does not effect 
equations for x„. 



Proof. The right- and left-hand sides of (|3.13|) are double-periodic functions of z that are 



holomorphic except at z = 0, where they have the form 0(z~^) exp((a;^_,_]^ — xD({z)). Such 
functions are equal if and only if the corresponding coefficients at z'"^ and are equal. 
The equality of the coefficients at z~'^ gives 



(3.17) 



which is fulfilled due to ( |0| ) (the second equality in ( p.l7| ) holds because f (x, t) is an elliptic 
function of x and, therefore, a sum of its residues is equal to zero). 

The equality of the coefficients at z~^ in the expansion of ( ^.13| ) at z = gives 
dt{\nh^^^,) - - x^)C«+i - xi) = - Mi+ 



+ ^ ^n+l Cix„J_^ X, 



-^n+lJ ^ '^n/'J Z^'^n S -^nJ ^ SK-^n -^n) ■ (3-18) 

The second line in ( |3.18|) is equal up to the sign to the sum of residues at x^, k ^ j, and at 
Xn+i, k ^ i, of the elliptic function 

Vnix, t) = Vnix, t) [C«+l ~ x) + ({x - X^) . 

Therefore, it equals to the sum of residues of this function at xl^_^_i and x^ We have 

t) + res^,^ty„(x, t) = {hi - /i;,+i)C«+i - O + M^- 



res. 



(3.19) 



Substitution of the right-hand side of the last equality into ( |3.18| ) imphes (after the shift 
n + 1 n) 



§ = 7n - 7n-l + E 2/^nC« " X^) " E [^n+lC( 



X'n ^n+l) + ^n-lC(2^n ^n-l) • (3.20) 
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From (|3.9| ) it follows, that ( |3.20| ) can be rewritten in the form 

T— ^— j- = 9a;G'^(a;^) + dtGl^{xD, 

where the function 



(3.21) 



= a. + In n.,.a^(x-x^) j' = ^« " ^'^"i' 

depends on t through the dependence on t of a;^ and a„, only. By chain rule we have 

d 



(3.22) 



dt 



Therefore, 



= (1 - ij,)a.GU<) + dt (g^«) 

From (|3.7] ) it follows that 

G'i(x^) = G^(x^) = G„(xO, 5.G;,(x;,) = -d^Glixl) = 9.G„(xJ, 

where 

/a(x - x„+i)cr(x + x„+i)a(x - x„_i)cr(x + x„_i)\ 

G„(a;) = a„ + In ^- — ■ . 

\ a\x + Xn) ) 

Therefore, equations ( |3.24| ) for i = 1, 2, have the form 

dt (G„(x„)) - (x„ - l)^^^ 



Xr> 



Xn 1 



dt {Gn{Xn)) - {in + 1))«9.G 



in + 1 

Equations ( |3.26|) and (|3.27|) are equivalent to the equations: 

Xn = {xl - l)d:rGn{Xn), dt {Gn{Xn)) = dt lu (x^ - 1). 



(3.23) 



(3.24) 



(3.25) 



(3.26) 
(3.27) 

(3.28) 



The first among them, coincides with equations ( |1 . 1| ) for Xn, and the second one (compare it 
with ( |2.14|) ) is equivalent to the definition of 7„ by (|3.16|) . Lemma and Theorem are proved. 



4 Direct problem. Spectral curves. 



In this section we consider periodic in n solutions of equations ( |1 . 1| ) . 
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Lemma 4.1 Let = a;„+jv(t) be a solution of Then 



n=l 



[xl - l)a^{2xr. 



is an integral of motion, I = const, and the monodromy matrix 



N-l 



T{t,z)= YlLn{t,z) (4.2) 



n=0 



satisfies the Lax equation 

dtT=[Mo,T]. (4.3) 

Proof. If Xn{t) is periodic in n, then the corresponding matrix functions Ln{t, z) and Mn{t, z) 
defined by ( |3.15D satisfy the relations 

Ln+N = Ln, Mn+N = Mn-dt{\nl). (4.4) 

Therefore, equation ( ^.14| ), dtLn = Mn+iLn — LnMn, imphes that 

dtT = -dti\nI)T+[Mo,T]. (4.5) 

Note, that if dtl = then ( |4.5| ) coincides with ( |4.3| ), and therefore, the second statement of 
the Lemma follows from the first one. 



Equation ( |4.5|) implies that the function 

Piz) = Iit)itiT{t,z)) (4.6) 

is time-independent. 

Matrix entries of L„ are double-periodic functions that are holomorphic on Fq except at 
z = 0. Therefore, (tr T) is also double-periodic and holomorphic on Fq outside z = 0. In 
order to prove that this function is meromorphic on Fq, it is enough to note that L^ has the 
form 

Lnit, z) = gn+iLng~\ (4.7) 

where 



From ( |3.2| ) it follows that in the neighborhood of 2; = 

L^ = izr'Ll + Li + Oiz), (4. 

where 



^ \ i -|- Xn+l i -t- Xn+1 I 



and 



^1 _ ^ / 1 — in+l \f —({Xn+l—Xn) ~-({Xn+l + Xn) \ (4 10) 

2\ 1 + 

Xn+l I \ CiXn+1 + Xn) (^{Xn+l Xn) I 
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Therefore, 

tr T = tr n Lr,{t, z)] = ^"^(1 + 0(^)). (4.11) 

\n=0 / 

The last equahty shows that (tr T) is a monic elhptic polynomial P^{z). Therefore, at z = 
we have P{z) = I(t)z~^{l + 0{z)). Hence, I(t) is an integral of ( |1.1D because P{z) does not 
depend on z. Lemma is proved. 

Due to ([4.3|) the spectral curve T defined by the characteristic equation 

R{w, z) = det {w - T{t, z)) = - (tr T)w + det T = (4.12) 
is time-independent. 

Lemma 4.2 The characteristic equation ([j-l^) has the form j[ j . 

Proof. We have already proved that (tr T) has the form ( |1.5| ). The relation $(x, z)^{—x, z) = 
p{z) — p{x), which is equivalent to the addition formula for the Weierstrass cr-function, implies 

det L„(t, z) = 2^2(x^_^i - 1) [p(x„+i - x„) - p(x„+i + x„)] . (4.13) 

Therefore, though Ln(t,z) depends on z, its determinant does not depend on z. Hence, 
(det T) is also independent. As it does not depend on t, we identify A^^ in (p.. 61) with 



N-l 

A^^ = det T(t, z) = 2-2^ n (^n - 1) (Pi^n - x„_i) - p(x„ + x„„i)) = 2-2^e^, (4.14) 

n=0 

where H is the Hamiltonian of system ( |1 . 1|) . Lemma is proved. 

For a generic point Q of the spectral curve F, i.e. for a pair {w, z) that satisfies ( [4. 12]) 
there exists a unique solution C„ = (c^(t, Q)) of the equations 

a+i(t, Q) = Ln{t, z)Cn{t, Q), dtCnit, Q) = M^t, z), (4.15) 

such that 

Cn+Nit,Q)=wCnit,Q), (4.16) 

and normalized by the condition 

cj(0, g)$(-Xo(0), z) + alio, Q)$(xo(0), z) = l. (4.17) 

Remark. Normalization ( [4.17| ) corresponds to a usual normalization \l/o(0,0,(5) = 1 of the 
solution ^„(x,t,g) of (pD defined by (U). 



Theorem 4.1 T/ie coordinates c\{t,Q) of the vector-valued function Cn{t,Q) are meromor- 
phic functions on F except at the preimages P± of z = 0. Their poles 7i, . . . ,7Ar+i do not 
depend on n and t. The projections zi^js) of these poles on Fq satisfy the constraint 

N+l 

E^(7.)=0. (4.18) 

s=l 
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In the neighborhoods of P± the coordinates ofCn{t,Q) have the form 

cl{t,Q) = z^-xl±it,z)e^^'+^-^'^>-\ (4.19) 

4(t,Q) = z^-xl±{t.z)e^^'-^-^'^>-\ (4.20) 
where x\±{t-,z) are regular functions of z: 

x\At^z) = zxl4t)+0{z'), Xn,-it,z) = xl_it) + zxli^it) + 0{z'), (4.21) 
such that the leading coefficients of their expansions have the form: 

W = c(t)(l - X.), Xl+W = c(t)(l + in), c(0) = 1, (4.22) 

Xl4t) = Sn{t), xl-{t) = -Suit), (4.23) 
where functions s„ satisfy the relation 

Sn+l = 2~^{xl^^ - 1) [piXn+1 - Xn) - <p{Xn+\ + X^)] Sn- (4.24) 

Proof. Vector-columns S^^ and 5*^^^ of the matrix-function 



n— 1 



So' = Si,, Sn{t, z)=l[ Lm{t, ^), n > 0, (4.25) 



m=0 



are holomorphic functions on Fq except at z = 0. They satisfy the equation Sf^_^i = LnS^K 
Therefore, the Bloch solution C„ of ( |4.15| ) has the form 

C,,{t,Q) = h^iQ)Sl'\t,z) + h2{Q)S^^\t,z), (4.26) 

where hi{Q), Q = {w,z) G F, are the coordinates of the normalized eigenvector of the 
monodromy matrix T{z), corresponding to the eigenvalue w. They are equal to 

hiiQ) = h,{Q) = -^^{w - T^\z)\ (4.27) 

where T^^{z) are entries of the monodromy matrix and the normalization constant r{Q) 
equals 

r(Q) = Ti2(^)$(-xo(0), z) + {w- T'\z)) <l>(a;o(0), z). (4.28) 

The function r{Q) has the pole of degree + 1 at P+ and the pole of degree N at P_. 
Therefore, it has 2N + 1 zeros. 

Let us show that of these zeros are situated over roots of the equation T^'^{z) = on 
one of the sheets of F. Indeed, if T^^(z) = 0, then eigenvalues w{z) of the monodromy matrix 
are equal to T^^{z) or T^'^{z). Therefore, r = at the points Q = {T^^{z),z). Equations 
( [4.27] ) imply that Cn has no poles at these points. The poles 7^ of Cnit, Q) on F outside the 
punctures P± are the other zeros of r{Q) and do not depend on n and t. Let us prove now 
that they satisfy ( p^ . 
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The function r*{z) = r{Q)r{Q'^) with a : Q ^ as permutation of sheets of F, is a 
well-defined function on To with the pole of degree IN + 1 at 2; = 0. As it was shown above 
it is divisible by T^^(z). Therefore, the ratio r*(z)/T^^(z) is an elliptic function with the 
pole of degree + 1 at 2; = and zeros at the points -z(7s). Divisors of zeros and poles of 
an elliptic function are equivalent. Therefore, ( ^.181 ) is proved. 

From ([4.7|) it follows that the vector-function Cn = Qn^Cn is a Bloch solution of the 
equation Cn+i = LnCn- Let us first consider the neighborhood of the puncture P+, which 
corresponds to the branch w = z~^{l + 0{z)) of the eigenvalue of the monodromy matrix. 

The vector-function with the coordinates given by (^4.22|) satisfies the equation 

That implies that in the neighborhood of P+ 



n+l 



[Xn, where 



is defined in (14. 81 



the vector-function Cn{t, Q) has the form stated in the theorem up to a time-dependent factor 
f+{t, z). Substitution of ( |4.19| , [4.20| ) into the equation dtCn = MnCn shows that dtf = 0{z). 
Therefore, the analytical properties of Cn near P+ are established. 

Now we are going to prove by induction that at P_ equalities ( [4.19| , [4. 201 ), and ( [4. 231 ) 
hold. For n = they are fulfilled by the normalization conditions. Let us prove first, that if 
and (|4^ ) hold for n' < n, then 



2 Kjf 



iCiXn+1 + Xn) - CiXn+1 - S„ + Xn~ + X' 



In,— 



0. 



Indeed, the equation C„+i = -L„C„ implies that Cn+i at P_ has the form 

(1 Xyi^i)Kn \ I ^n+l\ 



n+l 



^1 + Xn+l)Kn I 



Hence, 



C 



N 




n+l 



^2n-N-l ( ~ Xo)K,r, 
I (1 + Xo)Kr, 



+ Z 



(4.29) 



(4.30) 



(4.31) 



If Kn 7^ 0, then the last equality contradicts the monodromy property Cjy = wCq = O(z^). 
Therefore, k„ = 0, and then ( |4.3CI| ) shows that Cn+i has zero of order ?t, + 1 at P_. Therefore, 
a step of induction for equalities ( 4.19 , 4.20| ) is proved. The same arguments show that if 
( [4.23| ) does not hold then the vector Cn has zero of order {2n — N) which again contradicts 
the relation Cn = 0{z^). 

Equalities (|4l9| , lOOp near P„ are proved, possibly up to a time-dependent factor /_(t, z). 
Their substitution into the equation dtCn = M^Cn shows that dtf- = 0{z) and completes a 
proof of (gJ^, WM ), and (^) . 

Let Cn{z) be a matrix formed by the vectors C„(t, Qi{z)), corresponding to two different 
sheets Qi{z) = {wi{z),z) of F. This matrix is defined up to permutation of sheets. From 
( 4.19| - 0^ ) it follows that in the neighborhood of z = 



Cn{z) = 
Therefore, 











^1 Xn)c Sfi. 
^1 -|- Xn)C 



+ 0{z) 



^"e~*«(^) 



det Cn = —2cSnZ + 0(z 



(4.32) 
(4.33) 
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and from the definition of C„ we have Cn+i = LnCn- Hence, 

Sn+l = SndetLn, (4.34) 

which coincides with ([4.24|) . Theorem is proved. 

The correspondence which assigns to each solution = Xn+N(t) of ( p..l| ) a set of 

algebro-geometric data {F, D}, is a direct spectral transform. The following statement shows 
that the results of Section 2, can be seen as the inverse spectral transform. 

Corollary 4.1 The solution 

t, Q) = ciit, Q)^{x - Xn{t),t) + clit, Q)$(x + x„(t), t) (4.35) 

of equation ^3. d[) is equal to '^n{x,t,Q) = c{t)ipn{x,t,Q), where ipn{x,t,Q) is the Baker- 
Akhiezer function corresponding to T and the divisor D of the poles of Cn', the factor c{t) is 
defined in ( U-SSj ). 

All the solutions Xn{t) of / i7j|j have the form x„ = |(x^ — x^), where a;^(t) are roots of 
the equation 

e{nU + xiXt)V + Wt + Z) = 0. (4.36) 

Here 6{z) is the Riemann theta-function corresponding to T; vectors U, V, W are defined by 
(^■^ , \2. f^ ; vector Z corresponds to the divisor D via the Abel transform. 

As follows from the Theorem 4.1, the function \E'„ defined by ( |4.35| ) has the same analytical 
properties on F as the function c{t)'ipn- Therefore, they coincide. Equation ( [4.36|) immedi- 
ately follows from the formula ( ^.1U| ) for ipn- 



5 Action-angle variables 



Until now we have not used the Hamiltonian structure of equations (|1 . 1| ) . Moreover, a'priory 
it is not clear why a system that has arisen as a pole system of elliptic solutions of the 2D 
Toda lattice is Hamiltonian. The general algebro-geometric approach which allows to derive 
a Hamiltonian structure starting from the Lax representation was proposed and developed 
in 



n m H 



The main goal of this section is to construct action-angle variables for ( |1.1| ). First of all, 
let us summarize necessary results of the previous sections. A point (pri,x„) of the phase 
space M. of the system defines a matrix function Ln{z) with the help of the formulae 

Ui=T'h^^^^^,-xlz), (5.1) 

12 1 2 1 + e^" 

1 — eP" 

This function defines the spectral curve F (with the help of ( [4.12| )), and the divisor D of 
poles 7i, . . . , 7Ar+i of the Baker- Akhiezer function Cn{Q) = {cl^{Q), c^{Q)) 

Cn+i{Q) = Lr.{z)Cn{Q), Cn{Q) = wCo{Q) , Q = {w , z) G T , (5.3) 
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normalized by the condition 

c^(0, QM-Xo, z) + c2(Q)$(a;o, z) = 1. (5.4) 



The divisor D satisfies ( [4.18| ), i.e. defines a point of an odd part of the Jacobian J^^'iV) G 



J(r), which is defined fiber of the projection 

{71, . . . , 7jv+i} G J(r) ^ 2a;0+ = ^ z(7s) G Tq, (5.5) 



s=l 



corresponding to 0+ = 0. All the fibers are equivalent and can be identified with the Prym 
variety of V. Note, that a shift of 0+ corresponds to the shift x — x + a for the solution 
( |L8D of 2D Toda lattice. 



The correspondence 

(p„, x^)eM^ {r, D e J^^(r)} (5.6) 

is an isomorphism. The coefficients («», A) of equation (|1.6|) are integrals of the Hamiltonian 
system ( |1.1D . Equations ( |1 . 1| ) on a fiber over F of the map ( [5. 6] ) are linearized by the Abel 
transform (|27 



The main goal of this section is to construct the action variables that are canonically 
conjugated to the coordinates (pi, ... , (pN-i, 4>- '■ 

N+l 

4>k = ^fc(7s), (j)- = (j)N - 4>N+1, (5.7) 
s=l 

on the Prymmian J^^{T). Note, that 0+ = 07v + 4>n+i- 
Theorem 5.1 The transformation 

{Xn,Pn) I ^ (0l,...,0Ar_l,0_;/l,.../Ar) (5.8) 

where Ik are a-periods of the differential dS = ln{A~^w)dz: 

Ik= i \n{K'^w)dz, (5.9) 

is a canonical transformation, i.e. 

N N-1 

^ dpn Adxn= J2 i^^k A S(j)k) + SIn a 5(j)-. (5.10) 

n=l k=l 



Proof. First of all, following the approach proposed in |TB|, we define a symplectic structure 



on A4 in terms of the Lax operator and its eigenfunctions. After that we will calculate it in 
two different ways which immediately imply (|5.1CI|) . 

The external differential 6Ln{z) can be seen as an operator- valued one- form on Ai. 
Canonically normalized eigenfunction Cn{Q) of Ln{z) is the vector- valued function on Ai. 
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Hence, its differential is a vector- valued one-form. Let us define a two-form a; on by the 
formula ^ 

uj = - {lesp^Vt + resp^fi") , (5-11) 

where 

^ =< Cl^^{Q)5Ln{z) A 5Cn{Q) > dz. (5.12) 
Here and below < ■ > stands for the sum over a period of a periodic in n function, i.e. 



n=0 

C*{Q) is the dual Baker- Akhiezer function, which is defined as a co- vector (row- vector) 
solution of the equation 

C:^,{Q)L4z) = C:{Q), C%{Q) = w-'C*{Q), (5.13) 

normalized by the condition 

C;iQ)Co{Q) = 1. (5.14) 

The form u can be rewritten as 

u = ^reso Tr < {C-l^{z)6Ln{z) A 6Cn{z)) > dz, (5.15) 

where Cn{z) is a matrix with the columns Cn{Qj{z)), Qj{z) = {z,Wj) corresponding to 
different sheets of T. 

Note, that C*{Q) are rows of the matrix C~^{z). That implies that C*(Q) as a function 
on the spectral curve is: meromorphic outside the punctures; has poles at the branching 
points of the spectral curve, and zeros at the poles 7^ of C„((5). These analytical properties 
are used in the proof of the following lemma. 

Lemma 5.1 The two- form uj equals 

N+l 



ijj 



^ 5z(7s) A51n^/;(7,). (5.16) 



The meaning of the right-hand side of this formula is as follows. The spectral curve by 
definition arises with the meromorphic function wiff) and multi- valued holomorphic function 
z{ff). Their evaluations ^(7^), -z(7s) at the points 7^ define functions on M., and the wedge 
product of their external differentials is a two- form on M.. 

Proof. The differential Q, defined by ( |5.12| ) is a meromorphic differential on the spectral curve 
(the essential singularities of the factors cancel each other at the punctures). Therefore, the 
sum of its residues at the punctures is equal to the sum of other residues with negative sign. 
There are poles of two types. 

First of all, Q has poles at the poles 7^ of C^. Note, that 6Cn has the pole of the second 
order at 7^. Taking into account that C* has zero at 7^ we obtain 

les^^n =< C:^,5L^Cn > ASzi^s)- (5.17) 
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From (|5.3| ) and ( |5.13|) it follows that 



< C*_^_i6LnCn > — < 



(N-1 \ /n.-l \ 

n L^\5L,,il[Lm]Co>={C*^5TCo), (5.18) 
m=n+l / \m=0 / 



where T is the monodromy matrix. Using the standard formula for the variation of the 
eigenvalue of an operator: Sw = Cq6TCo, we obtain that 



res^^i7 = (5 In w (7s) A (52;(7s). 



(5.19) 



The second set of poles of $7 is a set of branching points qi of the cover. The pole of C* 
at qi cancels with the zero of the differential dz, dz{qi) = 0, considered as differential on 
r. The vector-function C„ is holomorphic at q^. If we take an expansion of C„ in the 
local coordinate {z — z{qi)Y^'^ (in general position when the branching point is simple), and 
consider its variation, we get that 



6Cn 



dz 



-5^(g,) + 0(l). 



Therefore, 5C„ has a simple pole at g^. In the similar way, we obtain 

dw 



Sw 



dz 



5z{qi). 



Equalities (|]20|) and (|]2T]) imply that 



reSo.f2 = reSo- 



< Cl_^-^6LndCn > 



At qi we have dLn{qi) = 0. Therefore, in the way similar to (|5.18|) , we get 

6wdz 



resg^fl = resg^ 



{C*^6TdCo) A 



dw 



(5.20) 



(5.21) 



(5.22) 



(5.23) 



Due to skew-symmetry of the wedge product we may replace 6T in ( |5.23|) by {6T — 6w). 
Then using identities C^{6T — 6w) = 5C*^{w — T) and {w — T)dCo = {dT — dw)CQ we obtain 

reSg^Q = -resg^ {6C%Co) A 6wdz = res,^ (Cj^(5Co) A 6wdz. (5.24) 

Note, that the dT does not contribute to the residue, because dT{qi) = 0. 

Expansions ( [4.19| , [4.20| ) near the punctures imply that 

resp± {C*^5Co) A 5wdz = 0. (5.25) 

Therefore, 



7V+1 N+1 

^reSg. {C^6Co) A 6wdz = — ^ res^^ {C%6Co) A 6wdz = ^ (51nw(7s) A Sz{'ys 

Qi S = l S = l 



The sum of ( p.lQ]) and ( |5.26| ) gives ( [5.16[ ), because 

N+l 

2ti; = — ^ res^^f2 — reSq,^. 
Our next goal is to prove the following statement. 



(5.26) 



(5.27) 
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Lemma 5.2 The symplectic form given by l \5.11\ ) coincides with the canonical symplectic 
structure 



N+l 
n=0 



(5.28) 



Proof. Using the gauge transformation ( [4 .71 











we obtain 



-reso Tr < C„l^{z)6Ln{z) A 6Cn{z) + A (5/„C„ 



(5.29) 



where 5/„ = dgnQn^- From ( [4.32| ), using the equahty 

(C(a;„+i + Xn) - C{Xn+l - Xn)) 5s„ + S{Xn!- + Xn,-) = SnS {({Xn+l - Xn) - ({Xn+l + Xn)) 

which follows from (|4.28|) , we obtain that the first term in (|5.29|) 

Ji = reso < Tr (C^li{z)6Ln{z) A 6Cn{z)^ > dz 

is equal to 

Jl =< 6hn+l A 6 (C(X„ - Xn+l) + C{Xn+l + Xn)) > ■ 

2s„+i 

Equation ([4.34| ) implies 

J _^ 25hn+l A 5Xn+l _ 25hn+l A / p{Xn - Xn+l) + p{Xn + J^n+l) \ ^ 

^ - 1 - 1 V^(a;n - a;„+i) - p(a;„ + x^+i)/ 

The second term in (|5.29|) equals 

J2 =< reso Tr (C^_|i5L„ A 5/„C„) > d^; =< reso Tr A > dz 



(5.30) 



From definition ( |5.1| ) of L„ by direct calculations we obtain that 



26hn+i A 5a;„ / p(x ) + p{Xn + Xn+l) , 
J2 — < — TT^ : I — : ^ : r I > • 



(5.31) 



hl+l - 1 V Pi^r^ ~ ^"+1) ~ Pi^n + Xn+l), 

At last, the third term in ( ^.29|) is equal to 

J3 = - < reso Tr (CnliSfn+i A 6LnCn) > dz =< reso Tr {{5Ln)K^ A 5fn+i) > dz , (5.32) 
because 



J4 =< reso Tr (c^l-^^6fn+i A LJfn) Cn > dz = 0. 
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(5.33) 



In order to prove ( p.33D , let us note that at z = 
Therefore, 









~X'n. 



(5.34) 



(5.35) 



J4 = <TesoTT{C-l,6f^^,ALJf^)Cn>p{z)dz 
= < reso Tr (C-_^//°+i A C„ > piz)dz. 

The last term in ( ^.35| ) is equal to the sum of residues at the punctures P± of the differential 

< Q+i^/n+i A Ljf^C^ > p{z)dz, 

which is holomorphic on F outside the punctures. Hence J4 = 0. 
From ( |5.32|) by direct calculations we obtain, that 



_ 25K+1 A 5Xn+l ^ 

^3 -< — p — n — > 



(5.36) 



Equations ( |5.30| , |5.3lD , and ( ^.361) imply (|5.28|) . Lemma is proved. 

Now we are ready to complete the proof of the Theorem. Equations ( [5. 161 ) and ( ^.141 ) 
imply that 



Indeed, we have 



N+l 

uj = —5a, = ^2 6\n{A~^w)dz. 
s=i ■'^+ 



N+l N+l 

6a (51nw(7s) A 52(7^) - iV51nA A 5 ^ ^(7^). 

s=l s=l 



(5.37) 



(5.38) 



The last term in ( |5.38| ) equals zero on the fibers 0+ = const of the map ( |5.5| ). 

The differential dS = \n{A~^w)dz is multi-valued on F but, following the arguments 
of one can show that its derivatives with respect to Ik, k = 1, . . . , N (which can be 
considered as coordinates on a space of curves given by ( |1.6|) ), are holomorphic differentials. 
dS is an odd differential with respect to the permutation of sheets of F. Therefore, In+i = 
—In and the definition of Ik implies that 



-^ds = dnl k 

dh 



Equations ( [5.7[ ) and 



l,...,iV-l, 
imply that 



■j^ds = dQ% - dn%^,. 

01 N 



(5.39) 



a 



N-l 

{(PkSik) + 

k=l 



-SI 



N, 



(5.40) 



and using ( [5.37| ) we finally obtain ( [5.101) . Theorem is proved. 
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